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Abstract. Using standard group theoretical techniques we construct the exact wave func- 
tion of the [70, 1~] multiplet in the orbital, spin and flavor space. This symmetric wave 
function is compared to that customarily used in the 1/Nc expansion, which is asymmetric. 
The comparison is made by analyzing the matrix elements of various operators entering the 
mass formula. These matrix elements are calculated by the help of isoscalar factors of the 
permutation group, specially derived for this purpose as a function of Nc. We also compare 
two distinct methods used in the study of the [70, 1~] multiplet. In the first method the 
generators are divided into two parts, one part acting on a subsystem of Nc — 1 quarks 
called core and another on the separated quark. In the second method the system is treated 
as a whole. We show that the latter is simpler and allows to clearly reveal the physically 
important operators in the mass formula. 
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1. INTRODUCTION 

For fifteen years the l/Nc expansion of QCD, wliere Nc is tfie number of colors 
has revealed itself to be an interesting and powerful approach for studying baryon 
spectroscopy. The method is based on an exact contracted SUc(2A''£) symmetry 
appearing in the large Nc limit, Nf being the number of flavors [sj, For large 
Nc this algebra becomes the S\J{2Nf) of the constituent quark model. Much work 
has been devoted to the ground state baryons where the operator reduction rules 
simplifies the expansion [6|. Usually higher order corrections of order 1/N^ are 
neglected. 

For excited states the problem is more complicated. To include orbital exci- 
tations, by analogy to the quark model, one can classify the large Nc baryons 
according to an extended symmetry given by the direct product SU(2iV/) xO(3). 
The group 0(3) implies the introduction of a spin-orbit and a tensor interaction. It 
is a phenomenological fact that these contributions are small so that the breaking 
of this symmetry is also small. An open problem is to investigate the validity of 
the l/Nc expansion in this extended symmetry. 

In the language of the quark model the excited states can be grouped into 
excitation bands A^. In the l/Nc expansion, the baryon masses have been calculated 



for the lowest multiplets of all excited bands from = 1 to 4. In these bands the 
multiplets belong either to the symmetric [56] or to the mixed symmetric [70] 
representation of SU(6). The symmetry of the wave function of excited baryons 
belonging to [56] representation allows a similar treatment as that of the ground 
state. The spin-flavor part being symmetric, the introduction of a symmetric orbital 
part does not modify the procedure. In the [70] representation (mainly the [70, 1~]), 
the situation turned out to be more complicated. There is a standard scheme 
[3] where the wave function is written as a product of a written on spa6.log. 
[stancu@ins symmetric ground state core composed of A^c ~ 1 quarks and an excited 
quark. In this approach, based on a Hartree picture, the s^^'^p orbital part is not 
properly symmetrized and the excited quark is always the last quark. The flavor- 
spin part is also asymmetric and corresponds to a single term of the exact wave 
function. 

Recently, a new scheme, which avoids the separation into a core and an excited 
quark has been suggested [i] . In that case the system is treated as a whole and the 
orbital-flavor-spin wave function is symmetric under any permutation of quarks. 
Some convincing quantitative arguments in favor of this procedure can be found 
in Ref. [9| The exact [70, 1~] wave function was written as a product of a core and 
a separated quark. Its orbital and the spin-flavor parts are mixed symmetric such 
as to recover the exact symmetric orbital-spin-flavor wave function. The procedure 
is described in Sec. 4.1. 

Using group theoretical arguments here we examine the relation between the 
exact and the customarily used asymmetric wave function. We argue that the 
description of the system is unsatisfactory when the spin operator S"^ and the 
isospin operators are separated into independent parts in terms of operators 
acting separately on the core and on the excited quark. Much better results are 
obtained when we directly consider the operators S*^ and acting on the whole 
system. In addition we examine the role of an operator constructed from the 
product of S, T and G generators of SU(4). 



2. SU(4) GENERATORS AS TENSOR OPERATORS 

The SU(4) generators Si, Ta and Gia, globally denoted by Eia are compo- 

nents of an irreducible tensor operator which transforms according to the adjoint 
representation [211] of dimension 15 of SU(4). We recall that the SU(4) algebra is 

[>S'i,TQ,] = 0, [Si,Gja] = i^ijkGkay 

[Tai Gil,] abcG ic) 

\Si,Sj\ = iSijkSk, [Ta,Tj,] =i£abcTc, 

[Gia,Gjb\ = jSij£abcTc~\~ jSab^ijkSk- (l) 

As one can see, the tensor operators Eia are of three types: Ei {i = 1,2,3) which 
form the subalgebra of SU(2)-spin, Ea (a = 1,2,3) which form the subalgebra of 
SU(2)-isospin and E^a which act both in the spin and the isospin spaces. They are 



related to Si, Ta and Gia {i = 1,2,3; a = 1,2,3) by 



E- = —- E = —■ E- 



(2) 



The matrix elements of every Eia between states belonging to the representation 
[Nc — 1,1] are given by 



{[N,-l,l]I'I'^S'S'^\E,am-l,l]IhSS^) = v/C[^-i.il(SU(4)) 





I'\([N^-l,l] [211] 
R ]\ SI S'l"' 




,(3) 



where Cl^=-i'il(SU(4)) = Nc{?,Nc + A)/S is the eigenvalue of the SU(4) Casimir 
operator for the representation [Nc — 1,1]. The three factors in the second line are 
respectively an SU(2)-spin Clebsch-Gordan coefficient (CG), an SU(2)-isospin CO 
and an SU(4) isoscalar factor. The necessary iscocalar factors for the derivation of 
the matrix elements of Eia have been calculated by Hecht and Pang p^] . Here, the 
phases and notations have been adapted to our problem. 

3. THE MASS OPERATOR 

The mass operator M is defined as a linear combination of independent operators 
Oi 

M = Y,CiOi, (4) 

i 

where the coefficients q are reduced matrix elements that encode the QCD dynam- 
ics and are determined from a fit to the existing data. Here we are concerned with 
nonstrange baryons only. The building blocks of the operators Oi are the SU(2A^j) 
generators Si, Ta and Gia and the S0(3) generators Their general form is 



(fc) ^(fc) 



n-l 



SFi 



(5) 



where O^'^^ is a fc-rank tensor in S0(3) and O^"^ a /c-rank tensor in SU(2)-spin, but 
invariant in SU(A^/). Thus Oi are rotational invariant. For the ground state one 
has = 0. The excited states also require k = l and k = 2 terms. 

The spin-fiavor (SF) operators O^gp are combinations of SU(2A^j) generators, 
the lower index i in the left hand side of (I5l) representing a specific combination. 
Each n-body operator is multiplied by an explicit factor of 1/N^~^ resulting from 
the power counting rules. Some compensating Nc factors may arise in the matrix 
elements when Oi contains a coherent operator such as G*" or T". 



3.1. The symmetric core plus excited quark procedure 



So far, for the baryons belonging to the [70, multiplet the general practice 
was to consider that they consist of one distinguishable excited quark moving in 
the collective potential generated by Nc — l ground state quarks [7], the latter 
subsystem being called core. The wave function of the core is symmetric both in 
the orbital (O) and flavor-spin (FS) spaces, which makes the treatment of the core 
analogous to that of ground state baryons. This description is known as the Hartree 
picture. 

To proceed, one defines separate SU(2A^/) generators that act on the excited 
quark s\ and and S^., and G*" that act on the core. Thus one has 

Si = s' + Sl T^ = e + T^, G'^ = g'^ + G';^. (6) 

As a consequence, the number of linearly independent operators Oi increases 
tremendously and the number of coefficients q, to be determined, becomes much 
larger than the experimental data available. For example, for the [70, 1~] multiplet 
with Nf = 2 one has 12 linearly independent operators up to order 1/Nc included 
0. For example, there is one operator of order A'^^: Net, three operators of order 
N^: i ■ s, 1/Nc i-t-Gc, 1/A'c i^^^ -g-Gc and 8 operators of order N'^: l/N^ t ■ T„ 
l/N, i ■ S,, i-g-T,, Si s ■ S,, £(2) .s-S,, l/Nl £(2) . t . {5„ G,} 
and C.- g- {S'c,Gc}- Then, in making the fit to the data, one faces the difficult 
situation of selecting among them the physically most dominant operators. We 
recall that there are only 7 nonstrange resonances belonging to this band. So one 
must select 7 out of 12 operators. Consequently, in selecting the operators one risks 
to make an arbitrary choice A much simpler method can be found, as shown 
below. 



3.2. A simpler procedure 

A simpler procedure is to avoid the splitting of the generators and the decoupling 
of the wave function and to consider instead only the global generators S**, T" and 
G*'* acting on the whole system of N^. quarks. However, the approach is not free of 
difficulties as the derivation of the matrix elements of the operators is more involved 
for a mixed symmetric wave function. Presently, the study of strange baryons is not 
possible. In the case of three flavors, one needs the analogue of Eq. ([3]) containing 
the corresponding SU(6) isoscalar. These factors have not been calculated yet. 

In addition to the fact that it uses an exact wave function, this approach implies 
only seven independent operators up to order 0{1/Nc) appearing in the mass 
operator: the order Ac operator three operators of order 1, £ ■ s, 1/Ac f'^^-' -G-G 
and l/A^c^-T-G and three operators of order C(l/A^c), namely l/N^ ^^ l/N^T'^ 
and 1/N^S-T- G. 



4. THE EXACT WAVE FUNCTION 



The total wave function is the product of the orbital (O), the spin (S), the flavor 
(F) and the color (C) parts. The color part being always antisymmetric, in order to 
fulfill the Fermi statistics, the orbital-spin-flavor must be symmetric. As the mass 
operator does not involve color operators, the color being integrated out, we are 
concerned with the orbital-spin-flavor part only. 

Here, as we are interested in the multiplet [70, 1^], the orbital and the spin-flavor 
parts must have both the mixed symmetry [iVc — 1,1]. In terms of inner products 
of the permutation group the wave function takes the form 

I mi) = E I [Nc - 1, l]Y)o\[Nc - 1, 1]Y)fs, (7) 

where Y is the corresponding Young tableau. Here we sum over the A'^c — 1 possible 
standard Young tableaux. The factor 1 / v^A^^c — 1 represents the CG coefficient of 
needed to construct a symmetric wave function [Nc] from its mixed symmetric 
parts. 



4.1. The decoupled wave function 

In order to calculate matrix elements of the operators listed in the Sec. 3.1 one 
must decouple the wave function as well. Using the Racah's factorization lemma, 
it is possible to decouple A^cth quark from the rest. The Sn^ CG coefficients can be 
factorized into an isoscalar factor times a CG coefficient of S^^^i. In the following, 
we need to know the position of the Acth quark inside a given Young tableau. In 
that purpose, one introduces the integer p which denotes the row where is the A'cth 
quark is located inside the Young tableau. 

The exact [70, 1~], but decoupled, wave function reads 



\iSJJs;Il3) = 

J2 a{p,icjq 



rrir. 



rria 



xK{[fW]p"\[N^-l,l]p)^^^^ 
X |4m,) l^.mi) |£,m,) |l/2m2) 1 1/2^2), 




(8) 



where ic and iq represent the angular momenta of the core and of the decoupled 
quark respectively and where a{p,ic,ig) are the one-body fractional parentage 
coefficients to decouple the A^cth quark from the rest in the orbital part. These 
are given by 



a(2. 







1) 



Nc 



(9) 



a(2,4 = l,^, = 0) = -y — , (10) 

a(l,4 = 1,^9 = 0) = 1. (11) 

The isoscalar factors K{[f']p'[f"]p"\[Nc — l,l]p) used in Eq. ([8]) are given in Ap- 
pendix A (Tables [10], [11] and [T2|) . The columns corresponding to p = 1 have been 
derived in Ref. [9|]. If we compare Eq. (JHl) with Eq. (3.4) of Ref. one can no- 
tice that in the latter only the terms with p = 2 have been taken into account. 
Furthemore, as the core was assumed to be in the ground state, the authors had 
considered a(2,4 = l,^g = 0) = and a(2,£c = 0,£g = 1) = 1. Thus the wave func- 
tion of Ref. [7| breaks Sn^ symmetry. As it represents only one part from the exact 
wave function we shall call it approximate or asymmetric. 

Tables [1] and [2] show the matrix elements for some spin and the isospin operators 
respectively calculated with the exact and with the approximate wave function. One 
can notice that the analytic expressions are different. Consequently, one expects 
the Cj coefficients determined from the fit to the data to be different if we use the 
exact or the approximate wave function. 

TABLE 1. Matrix elements of the spin operators calculated 
with the approximate (Ref. [3|) and the exact, Eq. ((Hj), wave 
functions, with s and Sc defined by Eq. 





Approx. w.f. Exact w.f. 


Approx. w.f. Exact w.f. 


48 
210 


N,+3 3(Afc-l) 

1 3(7V,-5) 

2 4Nc 

-, 3(Arc-l) 


Afc+3 3(^0-1) 
r, 3(37Ve-5) 

o 3(Afe-l) 
^ 27Ve 



TABLE 2. Matrix elements of the isospin operators calculated 
with the approximate (Ref. [3]) and the exact ,Eq. ([5]), wave 
functions, with t and Tc defined by Eq. ([5]). 







(T2) 




Approx. w.f. Exact w.f. 


Approx. w.f. Exact w.f. 



28 


Wc+3 


3(Arc-l) 


Wc+3 


3(^0-1) 








27Vc 


2Na 


48 


-1 




2 


3(Afe-l) 
2Na 


210 


1 


3{Na-5) 


2 


3(3iVc-5) 


2 


4Nc 




2Na 



4.2. The global wave function 

As already mentioned above, one can write the exact [70, 1~] states without 
decoupling them into a core and an excited quark. If there is no decoupling, there 
is no need to specify Y, the matrix elements being identical for all F's, due to 



Weyl's duality between a linear group and a symmetric group in a given tensor 
spac^. Then the explicit form of a wave function of total angular momentum 
J = i + S and isospin I is 



E ^ 

mi,ms 



i s 



J 

J3 



I [N, - l,l]ime) I [N, - 1, l]SmJh), (12) 



each term containing an SU(2) CG coefficient, an orbital part | [A^c — 1, Ij^m^) an a 
spin- flavor part | [Nc — 1,1] Smglh) . 



5. RESULTS 

Here we present the results obtained from different fits to the experimental 
data. In the fits, the seven nonstrange resonances have been taken into account: 
27Vi/2(1538±18), ^iVi/2(1660±20), ^N3/2{1523±8), ^Ar3/2(1700±50), ^A/'5/2(1678± 
8), 2Ai/2(1645±30) and 2A3/2(1720 ±50). 

In a first stage, we describe the fits obtained when we use the exact decoupled 
wave function. Afterwards, the results obtained with the global wave function are 
presented. 

In each case, we follow the spirit of the Hartree picture which leads to a one-body 
spin-orbit operator i-s. Its matrix elements are naturally of order N^. 



5.1. With the decoupled wave function 

Tables [3H6] show the four different fits considered. Each time, the results obtained 
with the exact decoupled wave function are compared to the ones obtained with 
the approximate wave function. 

In Table [3l we decouple the spin and the isospin operators. The operator 
1/NcTc-Tc is not present because its matrix elements are identical to of those 
1/A^c Sc ■ Sc for the approximate wave function (see Tables [H and [2]). This is 
apparently a practical advantage in the decoupling scheme but it has considerable 
physical disadvantages. One can notice that even if the Xdof is satisfactory, the fit is 
very bad. Indeed, the value of Ci is under-evaluated with respect to the commonly 
found value of around 500 MeV and the values of C3 and C4 are exceedingly large 
and of opposite signs, which suggest some compensation. 

The fits presented in Tables H] and [5] seem better. Here the linear combinations 
2s ■ 5c + 5c ■ 5c + 3/4 = 52 or 2t ■ + Te ■ + 3/4 = have been introduced. The 
coefficient ci has recovered its common value and the coefficients C3 or Cg have 
reasonable sizes, being about 70 MeV smaller for the exact wave function than for 
the approximate one. The quark-core operators 1/NcS-Sc or 1/Nct-Tc are still 



1 see Ref. [HI, Sec 4.5. 



problematic because the value of their respective coefficients are too high, of order 
500 MeV. 

The last fit shown in Table [6] correct these problems. All the coefficients have 
their natural sizes. This shows the necessity to consider the isospin-isospin operator 
on the same footing as the spin-spin operator. The values obtained with the exact 
wave function and the approximate one are identical in this case because the matrix 
elements of the operators considered are the same for the two wave functions. By 
construction, in both cases they are eigenfunctions of the total spin and isospin 
operators. 

TABLE 3. List of operators Oi and coefficients Ci obtained in the numerical 
fit to the 7 known experimental masses of the lowest negative parity resonances 
(see text). For the operators O3, O4 and O5 we use the matrix elements from 
Tables H] and H 



0, 




Ci(MeV) with approx. w.f. 


Q(MeV) with w.f. Eq. ^ 






211±23 


299 ±20 






3±15 


3±15 


03 = 


* '-'c 


-1486 ±141 


-1096 ±125 


04 = 


1 qi qi 
'-^c'-'c 


1182 ±74 


1545 ±122 


05 = 


1 j-arpa 


-1508 ±149 


417 ±79 


Xdof 




1.56 


1.56 



TABLE 4. Same as Table [3] but for O3, which combines O3 and O4 instead of using them 
separately. 







Ci(MeV) with approx. w.f. 


Cj(MeV) with w.f. Eq. ^ 


Oi 


= Ncl 


513±4 


519±5 


02 




3±15 


3±15 






219±19 


150 ±11 


O5 


1 farpa 


417±80 


417±80 



Xloi 1-04 1.04 



TABLE 5. Same as Table |4] but combining isospin operators instead of spin operators. 



0. 




Ci(MeV) with approx. w.f. 


c,(MeV) with w.f. Eq. jS]) 


Oi = 




516±3 


522 ±3 


02 = 




3±15 


3±15 


03 = 


J_ i qi 


450 ±33 


450 ±33 






214 ±28 


139 ±27 


Xdof 




1.04 


1.04 



TABLE 6. Fit with global operators proportional to the SU(2)-spin and SU(2)-isospin 
Casimir operators acting on the whole system (see text). 



Oi Ci(MeV) with approx. wi. Cj(MeV) with w.f. Eq. ^ 



Oi=Ncl 




484 ±4 


484 ±4 


O2 = ts' 




3±15 


3±15 






150±11 


150±11 




fr-r- + |) 


139 ±27 


139 ±27 


Xdof 




1.04 


1.04 



5.2. With the global wave function 

With the simphfied procedure described in Sec. 3.2 we can analyse the role 
of every of the seven independent operators introduced there. Table [7] shows six 
different fits to the experiment. The operators O5, Og and O7 are normalized to 
allow their coefficients Cj to have a natural size. As already emphasized, the Fits 
1-4 indicate that the coefficients of O3 and O4 have similar values. 

The partial contributions and the theoretical masses obtained from the Fits 1 
and 6 are presented in Tables [S] and M respectively. From Table [H] one can notice 
that the isospin-isospin operator in A masses plays a comparable role to the spin- 
spin operator in A^* resonances. This was impossible to observe in the symmetric 
core + excited quark procedure where the isospin-isospin operators were always 
ignored, for reasons explained above. From Table [9] one can see that the operator 
O7, never included before, is dominant in all resonances except "^Nj. This is a new 
finding, to be algebraically understood. 

The operators O5 and Oq do not seem to play an important role because, 
in addition to the fact that their coefficients are small and have an error bar 
comparable to their central values, their removal from the fit does not deteriorate 
it too badly. This justifies the previous choice presented in Section 5.1 where O5 
and Oq were neglected. Of course. Fit 4 is identical to the one shown in Table El 

Table [7| does not include a fit with O3, O4 and O7 together. In our calculations 
we found that the simultaneous presence of O3, O4 and O7 leads to a Xdof ~ 2. In 
this case the coefficients C3 and C4 become a bit higher, of the order 270 MeV and C7 
becomes negative suggesting a possible compensation with the contributions of O3 
and O4. This suggests that, by construction, O7 contains part of the contribution 
of the spin-spin and isospin-isospin interactions. As mentioned above, the role of 
O7 needs more investigation. 

6. CONCLUSIONS 



In principle, both the core + excited quark (Sec. 3.1) or the global (Sec. 3.2) 
procedures are legitimate as long as they are combined with adequately constructed 
wave functions. The core + excited quark procedure was the first attempt to study 
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the physically dominant operators in the mass formula. Our conclusion is that the 
spin operator 1/NcS-S is dominant in N* resonances, that the isospin operator 
1/NcT-T is equally important in A resonances and that 1/Nc S - T ■ G plays 
about the same dominant role both in A^* and A resonances except for ^iVj, the 
contribution to the mass being of the order of 200 MeV in all cases. More work 
should be done algebraically in order to understand the role of 1/NcS-T-G. 

Moreover, we found that all operators containing the 0(3) generators bring only 
small contributions to the mass, from 4 McV to 42 McV. This finding is consistent 
with the constituent quark model assumptions about the feebleness of the spin- 
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orbit and the smallness of the tensor interaction. 

We have shown that the separation core + quark procedure fails to emphasize the 
role of the isospin operator. This is due to the inherent structure of the asymmetric 
ground state core + excited quark wave function 0] which leads to equal matrix 
elements for S*^ and T^. Then the remaining part of the isospin interaction 1/Nct-Tc 
becomes exceedingly large if included in the fit (Table 3) and it is not surprising 
that in all previous studies (see Ref. [9] for a review) it has been totally ignored. 

In conclusion, the simple procedure we advocate here brings much more physical 
insight into the study of the nonstrange [70, 1~] baryons in the 1/Nc expansion. It 
is an urgent need to determine the isoscalar factors of SU(6) for mixed symmetric 
representations [A^c~l)l] in order to extend Eq. (3) to SU(6) and apply it to 



strange baryons. 



A. ISOSCALAR FACTORS 

Here we reproduce the isoscalar factors needed to construct the exact decoupled 
wave function (see Eq. ([8])). Detailed information can be found in Refs. 0, [HI, [TJ 



TABLE 10. Isoscalar factors K{[f']p'[f"]p"\[f]p) for 5* = / = 1/2, 
corresponding to ■^S when Nc = 3. The second cohimn gives resuhs for 
p = 1 and the third for p = 2. 

[/>'[/>" [A^c-l,l]l [A^c-l,l]2 



2(Ar^-2) Y 4Af, 



fJVo + l Ng-l l o F Af^ + l jVc-1 1 -I _1 / jVc-1 n 



L2'2J-^L2'2J^ 2V Af,-2 " 





-1 




-2 




-1 




-2 



TABLE 11. Isoscalar factors K{[f']p'[f"]p"\[f]p) for 5 = 3/2, / = 1/2, 
corresponding to ^8 when Nc = 3. The second column gives results for p ^ 1 
and the third for p = 2. 



[/>'[/'>" 



[iV,-l,l]l 



[A^c-l,l]2 



[Na±3 jVe-31 



2 


2 




JVc+3 


Na- 


3 


2 


2 




iVc+3 




3 


2 


2 




Arc+3 


Wc- 


3 



1[ 



2 


2 




Wc + l 




1 


2 


2 








1 


2 


2 




Wc + l 




1 



]1 ^ 



(jV,-l)(jVe + 3) 

JVc(Afc-2) 

5(jVc-l)(jV<;-3) 
2^0(^0-2) 

' (jVe-3)(jVe + 3) 

2V 2Afe(Wc-2) 



TABLE 12. Isoscalar factors K{[f']p'[f"]p"\[f]p) for 5 = 1/2, / = 3/2, 
corresponding to ^10 when Nc = 3. The second column gives results for p = 1 
and the third for p — 2. 



[/>'[/'>" 



fiV.- 1,111 



f Afc + l JVc-1 1 9 r jVc+3 iVc-3 1 9 1 / 5(jVc~l)(jVc-3) 
L2' 2J^L2' 2J^ 2Y 2Afe(A'c-2) 



2 


2 


Afc + l 


7Ve-l 


2 


2 


No + 1 




2 


2 


Ar„+1 





1[ 



iVc-31 



2[ 
1[ 



2 


2 




JV^+3 




3 


2 


2 




Afc+3 




3 


2 


2 




JV^+3 


Na- 


3 



1 / (JVc-l)(Afc+3) 

2 V A'o(JVc-2) 



1 / (jVc-3)(jVc+3) 
2Y 2Afe(iVc-2) 



]2 
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